1 






THE DETONATIC*' 
SOLID EXPLOS. /. 



BY 

CHESTER ARTHUR Z.'M 



. • . ' ’ ‘ V * 




Library 

U. S. Naval Postgraduate Schoo! 
Monterey, California 




i 



t 



- i - 



THE DETC»JATION OF 
SOLID EXPLOSIVES 



by 

Chester Arthur Zimmennan 



A THESIS 

Presented to the Graduate Faculty 
of Lehigh University 
in Candidacy for the Degree of 
Master of Science 



Lehigh University 

1952 



- ii 



This thesis is accepted and approved in partial fulfill- 
ment of the requirements for the degree of Master of Science. 




i 



•i*V 





- iii 



tabu; of contij^ts 

Page No. 

Introduction 1 

Initiation and Propagation of Detonation 3 

Shock Viaves 7 

Detonation Waves 12 

The Chapman-Jouguet Condition 12 

The Detonation Zone 19 

The Hydrodynamic Theory 22 

Gaseous Explosives 23 

Solid Explosives 25 

Variation of Parameters - Abel Equation of State 26 

Equation of State Based on Energy Considerations 29 

Halford-Wilson-Kistiakowsky Equation of State . , 31 

Calculation of the Detonation Velocity and Related 

Quantities 35 

Composition of the Detonation Products 39 

Results of the Calculations A6 

Appendices 

A. Relation of Detonation Velocity to c + u at Various 

Points on the Rankine-Hugoniot Curve 52 

B, Secondary Effects 55 

References 57 

General Bibliography 59 



- 1 - 



INTRODUCTION 



The purpose of this paper is to examine and correlate the 
various physical and chemiceil relations proposed to explain the phenom- 
ena of detonation, in particxilar, the detonation of solid e:q)losives. 

In order to assure a reasonably comprehensive treatise, it is desirable 
at the outset to devote some attention to the correction of certain 
vddespread misconceptions regarding explosives and their behavior. Two 
distinct types of explosion processes, deflagration and detonation, are 
evidenced by materials referred to as explosives. It has been common 
practice, especially in the early literature on the subject, to distin- 
guish between the two on the basis of chemical composition of the 
explosive and its intended application. Deflagration is invariably 
associated with "low” explosives or propellants auid the continual employ- 
ment of certain explosives such as nitrocellulose and black powder for 
their propulsive effect has led to the erroneous association of chemical 
compositions with the type of behavior. Similarly the detonation process 
is associated with primary and "high" explosives and likewise is con- 
sidered to be the sole type of behavior evidenced by other chemical 
substances. Unfortunately the above classification is not borne out by 
facts. On the contrary, experimental results indicate that all explosives 
are capable of exhibiting both modes of behavior. In view of this it is 
necessary to establish some criterion other than chemical composition to 
enable one to distinguish clearly between the two processes. 

The deflagration process, in its broadest sense, can be con- 
sidered a rapid burning in which the essential rate-determining factor 
is the speed of the chemical reaction. That is, the rate at which the 
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deflagration proceeds is related directly to the kinetics of the reaction 
taking place in the burning zone. 

The propagation of deflagration is closely associated vdth 
transport mechanisms; either the transport of heat or of atoms and free 
radicals acting as chain-carriers or a combination of the two. As long 
as mechanisms such as these control the propagation, there is established 
a limiting velocity beyond which the rate of deflagratic»i cannot pass. 
Extensive investigations have been carried out in connection with the 
process of deflagration in a gaseous medium and several excellent treat- 
ments are available (1,2). 

In the case of detonation, the characteristic feature distin- 
guishing the process is the appearance of a pressure wave propagating at 
a high velocity* The rate at which detonation proceeds is calculable by 
means of hydrodynamic and thennodynamic equations without the necessity 
of reaction-kinetic data. The application of the classical hydrodynamic- 
thermodynaaiic theory requires no assumptions regarding chemical mechanisms 
other than the postulate that the reactions occur at such a rate so as 
not to be a limiting factor. 

This theory, which is reviewed in detail, is based upon 
e3q)resslons for the continuity of mass, momentum and energy across a 
so-called detonation front and upon the assumption that the stable 
detonation velocity is the minimum one compatible with the conservation 
equations. When supplemented by ^ua equation of state and thetmal data 
for the products of reaction, the theory permits the calculation of 
detonation velocity, temperature and pressure, as well as other parameters 



of interest 
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INITIATION AND PROPAGATION OF DETONATION 

The currently accepted hydrodynamic detonation theory, which 
will be discussed in a later section, unfortunately provides no informa- 
tion concerning the mechauilsm by \diich a stable detonation wave arises. 

It is well established by experience that detonation can be brought about 
by a wide variety of conditions. Including mechanical shock, high 
temperature, friction, and adiabatic compression of interstitial gas. 

Only recently, however, have successful attempts been made towaixls the 
establishment of what might be considered a generalized theory of ini- 
tiation. The work of S. J. Jacobs and D. P. MacDougall (3) at the 
Explosives Research Laboratory of the N.D.R.C. involving the controlled 
deflagration of readily detonable substances and that of F. P. Bowden 
and his associates (4) on the Isolation of factors responsible for 
thermal ignition provided a basis upon which such a theory could be 
proposed. The approach adopted by both G. B. Kistlakowsky (5) and A. R. 
Ubbelohde (5) is that of a preliminary deflagraticai evolving discontinu- 
ously into a detonation. 

Consider a mass of explosive particles in which a small 
perlpherial deflagration has been started. As the deflagration spreads 
inward, pressure gradients develop due to the restricted flow of the 
gaseous products through the interstices between the particles. As the 
interior pressure rises, the speed of the deflagration increases, result- 
ing in a faster decomposition and a further increase in pressure. This 
constantly increasing deflagration may evolve into a stationary combus- 
tion in which the flow of heat and material away from the reaction z(xie 
compensates for that produced by the reactions, or the deflagration rate 



tit 





- 4 - 



may continue to increase until all the material is consumed provided the 
mass of explosive involved is less than the "critical size". The concept 
of "critical size" for a particular e^qalosive is a function of many 
variables including shape (volume-area ratio), particle size, crystal 
form for polymorphic substances, and density. It is probably best 
defined, although somewhat nebulously, as a limiting criterion, that is, 
the critical size is that size at which the rate of decomposition under- 
goes a discontinuity. The discontinuity may result in a partial or 
complete extinguishing of deflagration, with perhaps a re-initiation at 
an adjacent spot in the substance or it may be evidenced by the phenomenon 
of detonation. To see how the non-stationary deflagration might evolve 
into a detonation we must consider the flow of gases from a center of 
deflagration in more detail. As each crystal decomposes it acts as a 
center from which the gas molecules stream outwards, thus forming a 
pressure wavelet. If a nvtmber of these wavelets can integrate to give 
a streaming velocity corresponding to that accompanying a stable 
detonation wave for the system, we have conditions identical to those 
observed in a detonation. The type of flow just postulated has a very 
close analogy in the field of gaseous combustion. A column of gas acted 
upon by a piston subject to acceleration is frequently cited as an example 
of the mechanism of shock wave formation in gases. Further discussion of 
this mechanism is Included in the following section. 

Assuming a shock wave has developed it is of Interest to con- 
sider how its passage through the explosive can initiate the decomposi- 
tion. The details of the mechanism by which a shock with its discontin- 
uous rise in pressure and non-isentropic temperature effects activate 
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the explosive molecules is still the subject of considerable specula- 
tion. 

It has been proposed by Ubbelohde (5) that the most probable 
mechanism for the propagation of detonation is the collision of high 
speed atoms, radicals or molecules from behind the wave front with the 
quiescent layer of molecules just ahead of the front. Particles behind 
the detonation front are assumed to have an unusueilly high activation 
efficiency due to high random velocities at the detonation temperature 
and a high directed velocity resulting from the hydrodynamic conditions 
of propagation of a shock wave. Impact of a high speed particle is 
visualized as disrupting an unreacted molecule before the kinetic energy 
has had time to be distributed among all the vibrational degrees of 
freedom. If dissociation occurs after the redistribution of energy, 
then such a dissociation must be attributed to high temperature in the 
ordinary sense. Since redistribution requires a hhiaI T or zero activation 
energy, whereas dissociation requires an activation energy about equal 
to the strength of the bond broken, it is probable that redistribution 
of the kinetic energy will occur before dissociation in most cases. 

It has been suggested by Wendlandt (6), Eyring (7) and others 
that the mechanism of chemical activation depends primarily on the high 
temperature produced by the shock wave. Although considerable 
xmcertednty exists in regard to reactions at the extreme temperatures 
and pressures of detonation, it has been shown (8) that the temperature 
coefficient of the energy-releasing process is much smaller than the 
usual Arrhenius value. This, plus the fact that it is unlikely that 
sufficient time exists in the activation process for the vibrational 
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and rotational quanta to attain the full Maxwell-Bolt zman equilibrium, 
indicates that the temperature effect alone probably does not completely 
account for the activation. 

Still another theory has been proposed suggesting that the 
pressure in the wave front causes a bending of the molecules placing 
reactive groups in close proximity to each other, whereupon the 
decomposition occurs. It is felt that this theory is the least satis- 
factory of those presented since Bridgeman (9) has shown that high 
hydrostatic pressure alone has no effect on explosives. This theory is 
also untenable from the standpoint of chemical mechanism since practically 
every conceivable first step in the decomposition of organic nitro 
explosives is endothermic and the pressure activation concept requires 
an exothermic initial step. 

Thus at the present time it appears most reasonable to 
attribute the propagation of detonation to the energy possessed by the 
shock wave as a whole rather than to attempt a more detailed and perhaps 
erroneous description of the phenomenon. 
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SHOCK WAVES 



The theoretical basis for the development of a theory of 
detonation stems from theory of shock waves proposed by Hugoniot (10) 
in the nineteenth century. In mathematical language he formulated the 
difference between shock waves, or compressional waves of finite 
amplitude, and weak acoustic waves. The former were found to be 
pix>pagated with velocities which are greater than the acoustic velocity 
of the medium, the more so the higher the peak pressure in the crest of 
the wave. The medium itself is moving in the shocks, the wave velocity 
being slightly less than this mass velocity plus the acoustic velocity. 
Regardless of the original shape of the (finite) compressional wave in 
the medium, a shock front is formed eventually as the wave propagates 
through the medium. This means that pressure, density, and temperature 
suffer discontinuous increases in the front of the wave. In the course 
of its propagation through a medium a pure shock wave decays, this 
process being superimposed upon the normal weakening of waves due to the 
geometric factor of divergence, vinless mechanical energy is constantly 
fed into the shock wave. In other words, the mechanical energy of shock 
waves is being degraded and is then found as the heat energy of the 
medium after the passage of the wave. 

Since a detonation wave is basic ally a self-sustaining shock 
wave, a brief qualitative discussion of shock wave formation appears 
desirable. It is well established that in a compressional wave of 
finite amplitude the velocity of a portion of a compressional pulse is 
directly proportional to its amplitude. Those portions of the pulse 
with the greatest amplitude will travel faster than portions of lower 
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amplitude, there! ore changing the shape of the pulse as indicated in 
Figure 1, 




Formation of a Shock Front Pb>Pa 

Becker's treatment (11) of the origin of a shock wave in a 
gas is also of interest , Let a piston in a long tube be given a sudden 
small velocity. The gas in front of the piston will be compressed and 
a pressure wave will move ahead of the piston with the velocity of sound 
in the medium. The gas in front of the piston and behind the pressure 
wave will be moving with the velocity of the piston. If the velocity 
of the piston is increased by another small increment, a new pressure 
wave arises that moves somewhat faster than the first, since it moves in 
the gas already heated by the foregoing compression and moving with the 
velocity of the piston, Each succeeding wave will move more rapidly 
than the preceding one and will ultimately overtake the initial one to 
form a composite wave with an infinitely steep front; i.e., a discon- 
tinuity in pressure and temperature will be formed. 

The qualitative argument can be made more explicit by the 
Riemann formulation of the equations of hydrodynamics. Such a treatment 
is presented by Lamb (12) and the results confirm the conclusions reached 
above . 

By reversing the argument given for compression waves, it is 
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apparent tnat a rarefaction wave of finite amplitude, in which the later 
portions of the wave are regions of lower pressure, will spread out as 
it advances and a discontinuity will not arise. 

Returning to the shock front, which for the purposes of this 
discussion can be treated as a mathematical plane, it becomes desirable 
to establish equations relating the conditions on both sides of the 
discontinuity. The equations for the conservation of mass, momentum and 
total energy have been employed by Rankine (13) and Hugoniot (10) in the 
formulation of three conditions relating the pressui*e, specific volume 
and particle velocity benind a discontinuity moving with a velocity D to 
those properties in the undisturbed medium ahead of the discontinuity. 

Consider a shock front moving into a region in which the 
particle velocity is u^, the density eind the pressure from a 
region in which the same properties are denoted by U|^, ei and 



shock front 

Figure 2 

Conditions at a Shock Front Moving with Velocity D. 

In a time dt a mass of fluid ^Q(D-UQ)dt will enter unit cross-section 
area of the front and a mass ^2^(D-UQ^)dt will leave the front (see 
Figure 2). As dt is made very small so that the layers on either side 
are infinitesimal, the mass flow into the front must equal the mass flow 
away from it. Thus an e3q>ression for the conservation of mass may be 
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The change in momentum associated with the mass flow across 
the front must equal the net force per unit area. Thus for the con- 
servation of momentum: 

fo(D-Uo)(ui-u^) = Pi-Pq (2) 

Finally the conservation of energy requires that the net work 
doiteby the pressures Pj_ and Pq equals the increase in kinetic and 
potential energy. The work done by each pressure is Pudt for unit area, 
the kinetic energy is i u^ and the internal energy is E per unit mass. 
Therefore , 

Pl'^l”Po% = El-Eo+i('^i“'io) (3) 

The above equations may be rewritten in the form below by 
substitution of V = l/p and some algebraic manipulation. 

D = y (Pi-Po)/(Vo-Vi) (4) 

ui = y (Pl-Po)(VVl) (5) 

Ei-Eq = k (Po+Pi)(VVi) (6) 

The three equations above are often referred to as the Rankine 
Hugoniot conditions. 

Before proceeding to the consideration of detonation processes 
it is desirable to consider Equation 6 in more detail and to investigate 
the thermodynamic properties of the shock wave. Knowledge of the 
equation of state and the heat capacity enables one to calculate 
Ej^-Eq as a function" of P^ and so that Equation 6 becomes the Rankine- 
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Hugonlot equation relating the pressure and volume In a compresslonal 
wave of finite amplitude for given initial conditions. The relation 
describes the change in inteivial energy before and after compression by 
the shock wave as a function of the pressures and specific volumes. 

This change is adiabatic, since no heat is exchanged with the surroundings 
but it is not isentropic. The non-isentroplc nature of a shock (or 
detonation) wave was not appreciated by many of the early investigators 
and led to a conflict with the energy principle, A method of determin- 
ing the increase in entropy occurring across a shock front is discussed 
along with other detonation zone parameters in a later section. 
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DETQMATIQM WAVES 



The existence of detonation waves was originally reported by 
Berthelot and Vieille (14) in connection with studies of gaseous 
ejqplosions. Such a wave differs from a shock wave in two respects. In 
the first place, the wave is sustained and propagated through the 
mechanism of a chemical reaction supplying energy. Therefore, the 
chemical energy released must be Included in the calculations of the 
Rankine-Hugoniot conditions. The second difference is that the velocity 
of the wave is determined by the thermodynamic state of the material 
immediately behind the wave front rather than by the motion of an 
external boundary surface as is the case in pure shock waves. The three 
conditions relating properties of the material on either side of the 
detonation front in terms of the consez*vation of mass, momentum, and 
energy are no longer sufficient to determine the unknown quantities 
behind the front in addition to the velocity of the front. The problem 
differs from that considered in the previous section in that there is 
no controllable external condition, such as the velocity of a piston, 
available to determine one of the variables behind the front. 

The Chapman-Jouguet Condition 

An additional relation is required to make possible a unique 
solution of the hydrodynamic equations. Chapman (15) postulated a 
relation between the velocity of the detonation wave and the sonic 
velocity in the gas behind the wave. He also assumed that the minimum 
detonation velocity consistent with the Rankine-Hugoniot conditions pre- 
vailed in a steady detonation. Although his final results were in error 
due to assxaming that isentropic conditions prevailed at a discontinuity. 
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the basic concept was Ideutical to the noathematlcal egression published 
by Jouguet (16) two years later. The condition now bears the name 
Chapman- Jouguet and postulates 

D = Cj^+uj^ (7) 

where uj^ is the material velocity and cj^ the sonic velocity in the 
reaction products behind the wave front. 

Although a valid theoretical proof of the necessity of this 
condition is not yet available, there is little doubt that it correctly 
describes the obseirved phenomenon. Several plausible arguments have 
been advanced which aid in the acceptcince of this condition. The general 
form of the pressure-volume cui^e representing possible states of the 
products after detonation is shown in Figure 3* This curve is specified 
by the Rankin e-Hugoniot conditions, the thermochemical data concerning 
the explosive and the equation of state. 




The point A represents the initial conditions existing in the 
e^qplosive before the passage of the detonation wave. If the Initial 
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material velocity ia zero, as is invariably the case viith solid 
explosives. Equation 4 gives the detonation velocity in terms of the 
initial and final states. 

D = Vo ^ (Pi-Po)/Vo-Vi) (8) 

From Figure 3 it is evident that this relation can be written as 

D = Vo i tan 9 (9) 

where & is the an^e between the system of lines passing through point A 
and the (negative) V axis. For any value of 0 greater than 0]_, there 
will be two distinct final states (C and G) corresponding to the same 
Vcdue of D. For 0 - 0]^, the line AB is tangent to the curve at point B 
and is the final state corresponding to the minimum value of D. 

If the final state of the products corresponds to a point C 
on the R-H curve higgler than the point B, it can be shown (see Appendix A) 
that the velocity of the detonation wave is less than u+c and a reure- 
faction wave travelling at sonic velocity with respect to the products 
would move faster than the detonation and hence engulf and destroy It. 

(a rarefaction wave will arise from any process that abstracts energy 
from the reaction products such as heat loss or turbulence.) 

At the point B, the velocity of the detonation wave is equal 
to the velocity of sound in the reaction products plus the mass velocity 
of those products, thus a rarefaction wave proceeding with its maximum 
velocity will be unable to overtake the detonation wave. The mathe- 
matical proof that D = c+u at point B follows that given by Kistiakowsky 
and Wilson (17) • From Equations 4 and 5 the detonation velocity is given 
by 
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The velocity of sound in a gas is given by the relation (18) 

c -y (ap/ap)s « v(--pj^ 



( 10 ) 



(11) 



If there is a point on the Rankine-Hugoniot (R-H) curve 
having the same slope as the corresponding point on the isentropic curve 
(PV curve with constant entropy), then at that point 
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where the subscript r refers to the Rankine-Hugoniot curve. Differ- 
entiation of the Rankine-Hugoniot relation. Equation 6, holding Pq and 
Vq constant yields 



Now, from the first law of thermodynamics, dE * TdS-PdV, thus Equation 
12 can be written 
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It follows that the slope of the R-H curve through a point at which 
dSj^ * 0 can be ejqjressed as shown in Equation' 14 below 
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This is also the slope of the line AB which is tangent to the R-H curve 
at point B. Now the detonation velocity at any point on the R-H curve 
is given by Equation 4 and is 



D = V, 



' Pi-Pq \^ 

iVVl/ 



(4) 
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Therefore, at point B 

D = V. » is c, (»y Eq- li) (15) 

Combining this with Equation 10 yields 

D = u^fc^ 

The argument which is used to exclude points on the R-H curve 
below B is based on the entropy of the products. Consider the two points 
C and G lying on the same straight line from A and corresponding to the 
same detonation velocity D» The entropy of products at C is greater than 
at G since the transition from G to C is equivalent to passage from an 
initial state G across a shock wave of velocity D to the final state C, 
Since there is a definite increase in entropy accompanying the passage 
of a shock wave the entropy at C will be greater than at G, therefore 
the upper point will be a thermodynamically more probable state than G. 
This is true for all angles of 

Thus points above B are more probable than points below but 
points above B are not stable due to the effect of the rarefaction wave. 

The question of the validity of the Chapman-Jouguet condition 
has also been discussed by von Neumann (19) • By applying the conserva- 
tion Equations Uf 5 and 6 not only to the completely reacted material 
but also to each stage of the reaction, he has been able to substantiate 
the validity of the condition in certain cases. In addition to the 
single R-H curve in Figure 3 representing the final state of the 
products, there will be a series of R-H curves corresponding to the 
intermediate stages. Since the detonation velocity is the velocity with 
which the region of chemical reaction progress, it is the same for each 
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intermediate stage. Therefore, the line from the initial state P^, 
to each representing the condition of a given layer must have 

the same slope and must therefore be the same line for all degrees of 
completion of the reaction. Since in order to get from incomplete to 
complete reaction, the material must have traversed all intermediate 
stages, this line must touch all of the R-H curves corresponding to 
different degrees of completion. The family of R-H curves corresponding 
to the degree of completion of the chanical reaction from n = 0 to 
n = 1 has the general form shown in Figure 4. 




Figure 4 

Family of Rankine-Hugoniot Curves from n = 0 to n = 1. 

The motion of a point on the diagram representing the behavior of an 
element of explosive passing through the shock front and into a finite 
reaction zone must always remain on a straight line to satisfy the 
conservation equations. In the shock wave, the representative point 
must begin on the curve n = 0, at point A, must Intersect the curves 
for increasing n in order and finally end somewhere on the curve n = 1. 
In orxier that there exist a discontinuity in pressure and temperature 
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at the shock front, there must be another intersection of the straight 
line vdth the curve n = 0, corresponding to point C. Thus in the shock 
front the representative point runs from A to the limiting value at the 
point C where it again meets the curve n * 0. So far no chemical 
reactions have occurred but the explosive material has been “shocked” 
into a ccaidition where reactions can proceed. Now the representative 
point must proceed continuously from the curve n = 0 to n = 1; this 
requirement immediately restricts the possible intersections to points 
on n ■ 1 at or above B. If the straight line intersects the curve 
n * 1 at a point above B, the tangent to the curve at that point would 
be steeper than at B, so the velocity of sound in the products would 
exceed the detonation velocity, and a rarefaction wave would engulf and 
destroy the detonation wave. Thus the intersection with the curve 
n » 1 must be at B giving a point of tangency and corresponding to Chapman- 
Jouguet condition, D = u^-t-Cj^. 

Von Neumann edso discusses the case of intersecting R-H curves, 
but since this possibility does not correspond to any known physical sit- 
uation it will not be pursued further. 

By a straightforward thermodynamic procedure, Scorah (20) has 
sho)4n that when the Chapnan-Jouguet condition is assumed the ratio of 
entropy increase to available energy passing throiigh the shock front is 
a maximum. Although this result is in accord with the fundamental 
principle of maximum degradation of energy observed in natural processes, 
it is regarded only as a further substantiation of the condition rather 
than a rigorous proof. 
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THE DETONATION ZONE 



Before proceeding vdth the develojment of the fundamental 
theory, it is advantageous to have an idea of the order of magnitude of 
the dimensions of the zones in which the process of detonation occurs. 
The essential elements of a steady one-dimensional detonation process 
can be illustrated as in Figure 5* The detonation zone is considered 
to be stationary in space and the e^qplosive feeds into the zone from 
right to left at the detonation velocity D, 



^ Det, Zone > 





Reaction Zone 


Shock 

Zone 


Po> 
< 











Figure 5 

Steady Plane Detonation in a Solid Explosive 



Conditions throughout the zone are considered to be steady. The material 
first passes into a plane shock zone where abioipt Increases of pressure, 
temperature, mass velocity and entropy occur. This transition takes 
place within a narrow region whose length is determined by the viscosity 
and the thermal conductivity. If the dimensions of this shock zoie are 
very small in comparison with other lengths under consideration, the 
shock may be treated as a mathematical plane or discontinuity without 
error. For one-dimensional motion, the Navier-Stokes equations (21) of 
motion and energy for a steady state in which the shock wave is at rest 
are: 



u^ _ 
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and 

-dQ -i(D-u)2 _p (V -V) = k Xo ~ 

^ ^ D dx 

where ory^ is the viscosity, k is the thermal conductivity, u is the 
material velocity, T is the temperature and x is the cooi*dlnate in the 
direction perpendicular to the shock front. These differential equations 
have not as yet been solved in closed form for general values of the 
parameters. For the case of a pure shock wave ( -AQ ■ O) in liquid di- 
ethyl ether, Becker (22) has calculated the length of the zcne to be 
6,^A9 at a shock pressure of 10,000 atms. and 1.4A° at a pressure of 
100,000 atms. For a solid explosive with an assumed viscosity of 0.2 
centipose and a thermal conductivity of 10^ ergs cm“^ sec“^, I^yring (23) 
calculated the length of the zone to be of molecular dimensions. Calcu- 
lations including the effects of viscosity and thermal conductivity in 
the conservation equations have shown these effects to be of Importance 
only at the very front of the detonation wave, where the amount of chem- 
ical reaction is negligible. 

The distance between the points of Initiation of the chemical 
reaction and the point at which the Chapman-Jouguet condition is 
applied is called the reaction zone. There are a number of methods for 
estimating the length of this zone. Since each method rests on the 
interpretation of a different phenomenon, agreement among results found 
by different methods greatly enhances the accuracy of the determination. 
The results of (1) direct observation of the luminous zone, (2) extrapo- 
lation of low-temperature rate measurements, (3) surface reaction rate 
calculations for granular explosives, and (4) studies of the building-up 
of a low velocity wave to steady state conditions, all indicate that the 
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reaction zone for TNT under steady state conditions is about ^ centimeter 
in length. Although the exact figm*e vdll vary from e;qplosive to explo- 
sive, it is reasonably well established that the reaction zone length 
for solid ejq>losives is between 1 mm. and 1.5 cm. 

The time required for a particle to pass through the reaction 
zone can be calculated knowing the zone length, detonation velocity, and 
the material velocity. Reaction times of the order of magnitude of 
1 microsecond have been found for practically all organic e3q>losives. 
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THE HYDRODYNAMIC THEORY 



Application of the Rankine-Hugoniot relations and the 
Chapman* Jouguet condition to an actual chemically supported detonation 
pennits evaluation of the energy relation, Equation 6, in terms of heat 
of reaction and specific heat data for the reaction products. Since the 
shock zone is of molecular dimensions any contribution of heat by thermal 



be written 



conduction is neglected. The energy relation can now 

El-Eo • KPltP„)(Vo-Vi) = g ♦ K Vt.Up) 

"Tq oo '^1 

or introducing the mean specific heat at constant volume and assuming an 



dV (16) 



ideal gas 






(17) 



where Q is the heat absorbed in the reaction of M grams at the initial 
temperature T^. As indicated above the computation of conditions about 
the detonation front requires the knowledge of the following: 

1. Equation of state of the products 

2. Composition of the reaction products 

3. Heat capacity of the products 

4. Heat of reaction 

In the case of solid explosives, the selection of an equation of state 
for products at temperatures of the order of 3000 - 5000° A and at 
pressures of 50,000 to 200,000 atmospheres presents a formidable problem. 
Our present day knowledge of equilibrium constants vinder the above 
conditions is also quite meager. In view of this situation, it appears 
desirable to test the validity of theory under the less stringent con- 
ditions provided by gaseous explosions. 



- 23 - 



Gaseous .fclxplosives 

The approach adopted in this discussion is similar to that 
applied to mixtures of hydrogen and oxygen by Levds and Friauf (24). 
Since the equations used for ideal gases vdll later prove of value in 
the preliminary treatment of solid explosives, they will be fully 
developed here. The Rankine-Hugoniot conditions for conservation of 
mass and momentum. Equations 1 and 2, may be written as below, if 
Uq is assumed to be zero. 



D = 



Vo 




(18) 



pL r (19) 

’’O 

Elimination of from the above equations yields 

Pi-P. = p (V^i) (20) 

'>0 

At the ChafHnan-Jouguet point on the R-H curve the detonation velocity 
expression Equation 4 can be written 



D = 7 (Pi-Po)/(V„-Vi) = ^ (2X) 

For an ideal gas, the isentropic expansion law is 

Vi 

- constant 



Therefore, differentiating and substituting into Equation 21 

D = Vq/Vi ^nPiVi = Vq/Vi 7nniRTi/M (22) 

in which n^ is the number of moles of gas per M grams of product gases. 
Since the initial pressure Pq is negligible in comparison with Pq^ in 



Equation 20, substitution of Equation 20 into Equation 22 vdth the 



elimination of D yields 




(23) 



This e^qjression may now be substituted into the energy relation. 



Equation 17, again neglecting to give 




(24) 



If the chemical reaction under consideration went quantitatively and 
the distribution of atoms in the products was independent of temperature 
and pressure. Equation 24 could be solved for the final temperature T-j_ 
if the variation of and Cy with T were known. Actually, of course, 
both n^ and Q depend upon pressui*e and temperature. The calculation is 
therefore somewhat more tedious but the principle involved is unchanged. 

Cnee the final temperature is found, the detonation velocity 
can be computed from the combination of Equations 22 and 23 



Lewis and Friauf have calculated detonation velocities for mixtures of 
hydrogen and oxygen with various excesses of both components and with 
the addition of inert gas. They allowed for the eqxiilibria 

1 H 2 +O 2 - H 2 O, H 2 - 2H, and H20f^02 “ 20H. Although the dissociation of 

0 2 was not considered and no allowance was made for the excess heat 
capacity of oxygen molecules due to electronic excitation, there is no 
reason to suggest that a similar calculation undertaken today would 
differ significantly from their results. The following table shows 
their calculated values and the experimental results obtained by Dixon 



D = (Vi+1) J niRTi/riM 



(25) 
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(25) and Payman and Walls (26). 

Table I 



Mixture 

(2H2+0p) 

(2H2K)2) + 
(2H2*K)2) 
(2H2+O2) + 
(2H2+O2) + 
(2H2+O2) + 3 Np 
(2H2+O2) + 
(2H2+O2) + 
(2H2402) + 
(2H2+O2) + 



Detonation Velocity 
D, m/sec 





Ccxnputed 


Observed 


Deviation % 




2806 


2819 


—0.4 


102 


2302 


2314 


-0.7 


302 


1925 


1922 


0.2 


502 


1732 


1700 


2.0 


1N2 


2378 


2407 


-1.2 


3N2 


2033 


2055 


-1.1 


5N2 


1850 


1822 


1.1 


2H2 


3354 


3273 


2.5 


AHp 


3627 


3527 


2.8 


6H2 


3749 


3532 


6.1 



The agreement between the calculations and the observed measurements 
serves as an extremely good verification of the hydrodynamic-thermo- 
dynamic theory. With one exception the differences between the 
comparative values are all within the e:jq)erimental errors of velocity 
measurements. The detonation velocity of a nxunber of gaseous ejcplosives 
of the hydrocarbon-oxygen type also have been calculated and found to 
be in good agreement with experimentally determined values. 

Solid Explosives 

The extension of the fundamental theory to the detonation of 
solid e:!^losives requires the selection of an equation of state applicable 
to gases under extreme conditions of temperature and pressure. Use of 
the ideal gas law is immediately ruled out since the experimentally 
observed detonation velocities are quite sensitive to initial density 
and Equation 25 based on the ideal gas law shows no dependence on this 
factor. 



Another complicating factor arising with solid explosives is 
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the question of the composition of the reaction products. To consider 
all of the various equilibria possible when a typical explosive contain- 
ing atoms of carbon, hydrogen, nitrogen find oxygen detonates is a 
tremendous mathematical task and cannot be justified in view of the 
uncertainty of the thermodynamic functions of the possible products. 
Justification for ignoring the more complicated equilibria involving 
free atoms, polyatomic molecules and free carbon in approximate 
computations will be presented later. 

Variation of Parameters - Abel Equation of State 

To illustrate the form assumed by the hydrodynamic equations 
and to ascertain certain qualitative aspects of the process of detona- 
tion a particularly simple equation of state will be chosen as an 
example. An equation suitable for this purpose, but of less value in 
connection with practical calculations, is the Abel equation, 

P(V- o^) = RT, with the covolume oi assumed constant. 

The details of deriving the relations will for the most part 
be omitted. The derivations are quite straightforward and similar to 
those given by liyring and his associates (7). 

Evaluation of the sonic velocity appearing in the Chapman- 
Jouguet conditions. Equation 7, in terns of the equation of state yields 

(D-U3^)2 s (26) 

Vx“0< 

which when inserted into the Rankin e-Hugoniot relations and combinations 
of these relations gives the following five equations. (The initial 
pressure is negligible in comparison to P 2 ^ and is omitted.) 
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ir^ cx/Vf. 



(27) 



* 2R 



(28) 



Q-t-CyTo 2 y 



(29) 



= 



(Q+C^Tq) 2 (Y^ »1). 

(l+'>(^Vo)2 



(30) 



^ « Y ^Ac 



T+1 



(31) 



Equations 27 to 31 furnish the detonation velocity and all other 
properties of the reaction products provided we have data concerning 
the heat of reaction, the mean heat capacity of the products, the co- 
volume and the gas consteuit per gram. It is apparent that neither the 
detonation velocity nor any of the properties of the products depend 
upon the rate of chemical reaction. 

An estimate of the probable accuracy of the calculated 
properties can be made from the above equations. The final temperature 
is a function of the heat of reaction and the meaui heat capacity of the 
products and in particular is independent of the covolvuae; thus the 
temperature should be determinable to a higher accuracy than properties 
depending on additional data. In particular the final pressure is least 
accurately obtained since, in addition to dependence upon the thermal 
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data and the covoliune^ it is also a function of the gas constant which 
depends upon the mean molecular weight of the products. 

The properties of the reacting explosive can be determined 
at various intermediate stages by replacing Q in Equations 27 to 31 by 
nQ where n denotes the degree of completion of the reaction. It is 
interesting to note that such a procedure is another application of the 
mechanism used by von Neumann in his attempt to justify the selection of 
the Chapman- Jouguet condition based upon the family of Ranicine-Hugoniot 
curves. The variation of the properties, temperature, pressure, ratio 
of c+u to D and entropy of a typical explosive as the reaction proceeds 
is indicated in Figure 6. Each property is plotted in terms of percent 
of the final value. 




Figure 6 

Variation of Properties within the Detonation Zone. 
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At the shock front the pressure is tvd.ce the value attained at the end 
of the reaction zone, since the substitution of nQ into Equation 28 
and simplifications yields, to close approximation, the relation 

Pn » (1 + 

The mass velocity is alviays the same fraction of its final 
value as is the pressure therefore, the same curve represents the 
variation of both parameters through the reaction zone* 

The ratio of (c+u) to D exceeds 2 at the shock front and reaches 
unity at the end of the reaction zone, thus it is seen that v#ithin the 
reaction zone a rarefaction vrave of velocity c+u vdll continually over- 
take and decrease the velocity of a discontinviity progressing at a lov#er 
velocity. At the end of the reaction zone D = c+u, (the Chapman-Jouguet 
condition) and a stable detonation vrave is possible. 

The temperature and entropy shov^ sharp Increases at the shock 
front and a gradual rise throughout the reaction zone. 

Equation of State Based on Energy Considerations 

In order to improve the agreement betv/een calculated and 
obsearved values of the detonation velocity, equations of state more 
representative of the extreme conditions confronted have been suggested, 
Jones (27) has developed an equation suitable for use at high pressures 
based upon Bridgman's (28) results vrith nitrogen at pressures up to 
15,000 atmospheres. The equation takes the form 




vfhere Eo(V/N) is assumed to be the potential energy of interaction of 




I 
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molecules in the gaseous product and N is the number of moles. The 
internal energy is expressed in the form 

Eo(V/N) • Ae~ + B(V/N) + C 

in which the first term represents the repulsive forces and the second 
the attractive forces. The constant C is assigned a value to make the 
minimum energy equal to the heat of vaporization. The function f is 
evaluated from compressibility data and the constants A, B and •< from 
the energy of the normal modes of vibration of molecules in the fully 
excited solid state. 

At lower temperatures j Jones assumes an equation of the virial 
form, the constants being adjusted so that the pressure computed at cin 
intermediate value of V by the both relations is the same. Extending 
the interaction effects determined by Bridgeman to higher temperatiires 
and densities the total energy of the reaction products can be expressed 
as 

E = Z[NiEi + (Z%-Na)(Eo^-3/2RT)] 

where is the number of moles of each molecular species, both gas and 
solid, E^ is the energy per mole exclusive of the interaction of 
molecules in the gas phase and N. is the number of moles of solid 
products. The last term, therefore, accounts for the energy of molecular 
interaction plus the vibration energy of 3/2RT per mole. The energies 
Ej^ can be evaluated from specific heat data at ordinary pressures since 
the interaction effects at high pressures are represented by E^. To 
determine the composition of the reaction products as a function of 
temperature, the activities of the various products are introduced. 
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The activity of each species can be expressed in terms of the composition 
and the equilibrium constants for the assumed reactions as a function 
of temperature at normal pressures to which ideal gas conditions are 
applicable. The increase in internal energy of the products must be the 
sum of the chemical energy Q released by the reaction and the work done 
by the pressure at the detonation front. Thus a knowledge of E, Q and 
the equation of state provides sufficient information to solve the 
detonation front equations. 

Although the above procedure involves several major assxmptions 
to make it feasible, the calculation is significant because the results 
tabulated in Table III give an idea of the success obtainable using an 
equation of state fitted to data other than that obtained by detonation 
measurements . 

The Halford-Wilson-Klstiakowsky Equation of State 

Information on gas states at the pressures and temperatures 
encountered in detonations can be obtained by the use of the hydro- 
dynamic theory if the variation of detonation velocity is knovn as a 
function of the initial specific volume (loading density). Since the 
detonation velocities of solid eaqplosives show a definite increase with 
loading density, such measurements provide a suitable means for finding 
the deviations of explosion gases from ideal gas behavior. The detonation 
velocity-loading density relationship is of little value in suggesting 
the correct form of an equation of state, but, once a form is established, 
the relationship is sufficient for the quantitative evaluation of its 
empiricsil parameters. Unfortunately, the pressure-volume relationship 
is rather insensitive to the form chosen which suggests the measurement 
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of detonation temperature for determining the appropriateness of the 
choice. Since detonation temperatures have not as yet been determined 
to a degree of accuracy sufficient to provide a check, we are forced to 
accept the form which fits the data best in the largest number of cases. 

An equation of state found to be very successful in the cal- 
culation of detonation velocities is one developed by R. S. Halford in 
cooperation with Wilson and Kistiakowsky (17)* The equation of state 
has the following form when all the products of reaction are gaseous 

PVM = nRT(UxeP^) (32) 

with X * 2niki/T"^VM 

V is ^he specific volume, n is the number of moles of gas per M gramsj 
P, T and R are the pressure, temperature and the gas constant per gram 
mole, lespectivelyi and n^^ and kj^ are the number of moles and an empirical 
"covolume" parameter of the i-th species, respectively; the summation 
extends over all gas species present. Satisfactory results have been 
obtained with the constants ois o,25 and ^ = 0*30. 

In cases in which solid products are formed in the detonation, 
the equation is slightly modified. It is assumed that the volumes of 
gas and solid are aduitive and that the gas obeys Equation 32 independent 
of the presence of solid components. Since thermal expansion 5ind baric 
compression of solids have opposing effects in respect to the volume, and 
since these effects are hot known at the elevated temperatures and 
pressures encountered, it is assumed that the specific volume of a solid 
component is independent of temperature and pressure. Thus the specific 
volume of solid components may be assigned values determined under normal 
conditions or may be treated as adjustable parameters. With these 
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assumptions, the equation of state for a mixture of solid and gaseous 
components becomes 

PM(V-Vg) = nRTF(x) (33) 



where Vg is the volume of solid constituents per gram of mixture, 
obtained by 

Vs = ^a Va + ^b Vb + (33a) 



^a* '^b**** weight fractions and V^, Vb,..» are the specific 

volumes of solid constituents a, b,.,. F(x) is a function defined by 

F(x) s 1 ♦ xe (33b) 



with 

X = k/T*’ M(V-Vg) (33c) 

k = Zn^kj^ (33d) 



The value of ec and ^ and all other symbols are the same as those used 
with Equation 32. A theoretical evaluation of the covolume factors 
(c.f. Equation 33d) is impractical in view of the extremely high density 
encountered in the reaction zone. These factors are therefore determined 
by comparison of the results of a calculation of detonation velocity 
based on an assumed value with the actual observed velocity. As an 
initial approximation only the most abxmdant detonation products are 
considered and the total covolume constant k of the mixture is obtained 
for several explosives at a number of loading densities. To evaluate 
the specific covolume factors kj^, it is assumed that k is an additive 
function of specific values. The individual values are then adjusted 
by an iterative procedure for optimum conformance with the mixture values. 
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since the original work was completed, other methods of evaluating these 
constants have been suggested and a wider range of explosives has been 
used as a basis, yielding slightly modified values. It has been found 
that the calculation of product composition is much more sensitive to 
the individual values of than is detonation velocity. 

The application of the Halford-Wilson-Kistiakowsky equation 
of state to the calculation of detonation velocities and related 
quantities will be illustrated below. 




<jx,i 
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CALGULATIOM OF THE DETONATION VELOCITY AND RELATED QUANTITIES 



The calculation of the detonation velocity and related 
quantities using the HalXoi^l-V^ilson-Kistiakowsky equation of state takes 
into account the effects of gas imperfections and results in an entirely 
acceptable set of parameters in the vast majority of explosives to vfhich 
it has been applied. 

Calculation of Detonation Velocity 

From pure thermodynamics, we obtain the relation 




Differentiating the equation of state, Equation 33> holding T 
constant gives the result 




= - P(U2xe 4 (9 



(VJ 






(35) 



from Equation 33b, the coefficient of (V-Vg) can be represented 

by F(x) and that of P by y(x) where 

y(x) = F(x) f xdF(x) (36) 

dx 



Thus Equation 35 can be written 

Uvit (v-Vg) F(x) V\dV /3 



(37) 



Substitution into the basic equation for D, Equation 21, yields the 
expression 



D = y(x^)/ (V^-Vg) F(xi) 



(38) 
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Introducing the equation of state once again gives 



D 



^1-^s 






( 39 ) 



Equation 38 niay be substituted into Equation U to obtain an ejqpression 
in terms of the specific volumes (Pq is neglected) 



p - Piy(3Ci) yi(Vo“Vi) 
F{x^) (V;^-Vg) 



W 



which on rearrangement becomes 

Vi-Vg y(x^)y ^ 



(41) 



Calculation of the Specific Heat Ratio 

From thermodynamics, one obtains the following equations for 
the specific heat at constant volume and the ratio of the specific 
heats at constant pressure to that at constant volume 

Cy = Cj + (d^p/dT^)^ dV (42) 

where Cy is the ideal heat capacity per M grams and 



Y= 1 



2M /l^\ 

Cy \dT/y / Idv)^ 



(43) 



Introducing the equation of state in these expressions yields 
Cy = Cy + nRo<(z(x) - 1) 

and 'y = 1 + nR [z(x)VCyy(x)j 



(44) 

(45) 



where 



z(x) = F(x) -OCX 

dx 



(46) 
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Calculation of the Temperature 

For an imperfect gas, Equation 16 ejq>resses the application of 
the first law of thermodynamics 

M(Ei - Eo) = (Ti-To)+M dV (47) 

(The reaction is first carried out to the final temperature at a large 
volume, such that the product gases are ideal, and these gases are then 
compressed to their final volume with a resultant energy term due to the 
gas imperfections.) 

Introducing the equation of state and standard thermodynamic 
relations, one obtains 

M ^'^^”'^®^dE/dV)x, dV * n^RT^^{F{x^)-l) (48) 

Equation 41 can be rearranged to the form 

V„-Vi . (49) 

which when combined with the Hugoniot Equation 6 yields 

El-^o “ ^Pif(xi)Vi/y(xi)yi = in3^RTiF^(xi)/y(xi)'3'^M (50) 

Therefore, from Equations 47, 48 and 50, one obtains 

Ti = ^ To-Q/cJ«ilRo<(F(xi)-l)-nj^RF2(xi)/2y(xj^)Ti (51) 

Solution of the above equations is based on an assumed knowledge of the 
final composition, thus an iterative approach appears desirable. If a 
large number of explosives are to be investigated a two-stage method is 
recommended. This involves the calculation of ideal values of D, T, and 
Y based on the assumption that the product gases are ideal and then 




35' 
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evaluating the ratio of actual to ideal values using tables prepared 
for this purpose. 

The entropy increase occurring in the detonation zone can be 
estimated if values for the increase in temperature and volume across 
the shock wave are assumed. 

The entropy increase across the shock zone due to an irreversible 
compression to a pressure P and temperature T is given by 

The entropy change due to the irreversible chemical reaction 
can be calculated by proceeding from the initial to the final state by 
hypothetically reversible steps, applying the relation 

+ nCpln ^ - n R In ^ (53) 

io 

in which the subscript 1 refers to the conditions prevailing at the end 



of the reaction zone 
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COMPOSITION OF THfc DETONATION PRODUCTS 

In order to carry out calculations of the detonation velocity 
and other parameters of interest, it is necessary either to calculate 
the equilibrixun composition of the reaction products or to make an 
assumption regarding this composition. Before discussing the detailed 
methods used to detennine the composition, a few general remarks concern- 
ing equilibria among possible products appear desirable. 

Dissociation into free atoms and radicals is fortunately not 
of importance at the high pressures encountered in detonation of solid 
eaqilosions. The dissociation constant of H 2 = 2H at 5000°K has been 
calculated statistically to be 44.7 atm. Since the concentration of 
hydrogen gas rarely exceeds 10^ by volume in most e^qplosives, it is found 
that at a pressure 10^ atmospheres about of the hydrogen is 
dissociated into atoms. This corresponds to an absorption of heat of 
about 3 kilocalories per kilogram of explosive, whose total heat of 
explosion is of the order of 1000 kcal. per kg. Thus, the final 
temperature is lowered to the extent of 0,3% by the inclusion of this 
dissociation but at the same time, the number of moles of products and 
the pressure is increased. The two effects oppose each other as far as 
detonation velocity is concerned and it is estimated that the net effect 
results in a 0.10^ decrease in the velocity. The dissociations of O 2 
and N 2 occur to a much lesser extent and can also be ignored in all 
ceilculations. Although the dissociations of water: HgO « ^ H 2 + OH 

and H 2 O = ^ O 2 + H 2 take place to only a very small extent these 
equilibria are considered in precise calculations. 

Various polyatomic molecules have been detected in gas samples 



I 
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analyzed after closed b<»nb detonations. Schmidt (29) has reported small 
but variable amounts of CH^, C 2 H 2 > HCN and C 2 N 2 from almost all 

explosives. It is suspected that the majority of these molecules are 
formed after the detonation is complete and while the products are 
cooling. Allowance is made in precise computations for the formation of 
methane and ammonia, but the remaining three au:e neglected in the 
Interest of mathematical simplification. If it is assiuned that the 
analytical data represents the tz*ue conditions existing in the detonation 
wave, the formation of all these complex molecules has a small but not 
negligible effect on the results of calculation. For example, in TNT 
where the oxygen deficiency is large and the polyatomic molecules listed 
above are particularly abundant neglecting their formation results in an 
overestimation of the detonation velocity by about 

The formation of free hydrogen in absence of solid carbon is 
considered in precise computations even though the equilibrium constant 
of the "water-gas" reaction, H 2 O + CO = C02+‘H2 ranges from 0.2 at 3000° 
to 0.07 at 5000®K. Since the number of moles of products does not change 
the only effect on the calculations is that due to the evolution of heat. 
The neglect of this equilibria in the absence of solid carbon can result 
in an underestimation of the rate by only a few percent. If, on the 
other hand, free carbon is present, the reaction C+H 2 O = CO+H 2 proceeds 
almost quantitatively. The equilibrium constant of this reaction is of 
the order of 10^ atm. at 3000® and 10^ atm. at 5000®K. Since the 
fugacity of CO in the reaction zone is usually less than 10^ atm., the 
ratio ^H 2 /^H 20 greater than unity as long as carbon is present. 

The assumption that the reaction goes completely to the right results in 
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an overestimation of the detonation rate since the heat absorption is 
overcompensated for by an increase in the number of moles of gaseous 
products. 

As a first approximation results may be obtained using the 
follovdng arbitrary decomposition equations for e 3 q>losives of the 
general formula, CqHj.OgN^. 

Case A; (q r/2 4 S) 

CqHj.OgN^ = r/2 H 2 O t (s - q - r/2)C02 + (2q - s + r/2)C0 + t/2 N 2 
Case B: (q + r/2 S) 

C^H^OgN^ = r/2 H 2 O + (s - r/2)C0 (q - s + r/2)C + t/2 N 2 

These equations express the assumption that all the hydrogen 
is converted to water, that the remaining oxygen then converts carbon 
to carbon monoxide and that additional oxygen, if any, converts carbon 
monoxide to carbon dioxide. The results of a calculation using the 
arbitrary product composition given by the above equations are shown in 
Table III, 

If one assumes that the detonation products are in chemical 
equilibrium at the final temperature pressure and density, then the 
equilibrium relations together with the stoichiometric restrictions 
provide the necessary relations for the calculation of the product 
composition. The thermodynamic equilibrium constants for the significant 
equilibria can be computed from the function - (F®-E®)/T for each of 
the components, together with for the equilibria of interest. 

Tables of these quantities for H 2 , CO, CO 2 , H 2 O, N 2 , OH, 02» and NO 
eind for the following equilibria over the temperature interval, 2000® - 
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5000®K. have been published (30). 

H2 + CO2 = CO + H2O 
2C0 = CO2 + C 
H2O = iH2 + OH 
2H2O = 2H2 + O2 
H2O + ^N 2 = no + H2 

The problem is conveniently approached by the consideration of 
the work content A of M grams of a mixture of gases at a large volume 
V® and temperature From the thermodynamic relations, we have 

A(V®, Ti) = rniy«4 + RTi ln7T(pi)“^ - nRTi (54) 

o 

where nj|^ is the number of moles of the i-th constituent and y** ^ and 
pj are its chemical potential and partial pressure, respectively, in 
the standard state (pure i-th component at and unit pressure), 
n is the total nvunber of moles. 

Equation 54 can be written 

A(V°, T]^) = Xniyuf + [n In p/n 4 Irii In nj^ - n](55) 

The work content at a volume is obtained by 

A(Vi, T^) = A(V°, Ti) - M j^JpdV (56) 

To evaluate the integral requires the introduction of an equation of 
state. Since calculations reported in a later section have been based 
on the Halford-Wilson-Kistiakowsky equation. Equation 33 » this will be 
used. To simplify the resulting e3q)ression it is convenient to use the 
function F(x) as defined in Equation 33b, . Performing the indicated 
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integration and simplifying one obtains 

A(ViTi) = + RT]l {“i In + n4n [ Pj^/nFCxj^)] 

n(e^^ - 1) /p -n} (57) 



The chemical potential is given by 

-^ = ( 3 a / 3H\v^a <58) 

and therefore 

“ ^ In + ln(P 2 ^/n) 4 2,303 /H. (xj^) 4 

2.303 n ki / 2 (x)/k] (59) 

where fi (x) ■ 0.4343 (e - 1)/^ - log F(x) (60) 

/I (x) = 0.4343 k e?* (6l) 



These functions have been tabulated with x as the argument (31) • The 
chemical reaction may be represented by 

ZWji Xj^ = 0 (62) 



where the mole number Wj^ for the i-th constituent in the J-th reaction 
is taken positive for products and negative for reactants. Thus 
Equation 59 may be written 

^ji/'i ■ '^ji^i'**' ^"^l [^ji In nj^ 4 (In (Pj^/n) 4 

2.303 n.(xi) 4 2.303 n ZWj^ ki r 2 (x)/k] ( 63 ) 



Now, since at equilibrium = 0 and by definition = 



- RT In Kp^ we may write 



log - log Gj 



(6A) 



where Kp^ is the thermodynamic equilibrium constant for the j-th reaction 
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and the q\xantity Gj is defined by the equation 

-log Gj = [log(P3^/n) + Q fa 

ZWji k^/k (65) 

Thus one can now write for the equilibrium constant in terms of 
moles (Kj) 

Kj = Kp . Gj « 7T (66) 

J 

If there are 1 constituents in the mixture, composed of m different 
atomic species, there will be (1 - m) relations of the form of 
Equation 66 and m stoichiometric relations, A convenient iterative 
procedure and the tabulated functions necessary for the calculation of 
the equilibrium composition of the reaction products are available in 
reference (31)* The equilibria usually considered include the five 
equations given above as well as 

3H2 CO = CH4 + H2O 
and 3/2 H 2 + ^2 " 

A comparison of the final composition of some typical explosives 
based upon the arbitrary decomposition equation and the equilibrium 
considerations is given in Table II below. The composition is given 
as the number of moles of each constituent per mole of explosive. 
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Table II 





TNTfJ, 


5 = 1.0 


Tetryl - 1.0 


Cyclonite ^ = 1.3 




Arbit. 


Equil. 


Arbit. 


Equil. 


Arbit. 


iquil 


H2 


0 


0.09 


0 


0.08 


0 


0.02 


CO 2 


0 


0.18 


0 


0.24 


0 


0.16 


CO 


3.5 


3.38 


5.5 


5.27 


3.0 


2.78 




2.5 


2.26 


2.5 


2.24 


3.0 


2.87 


N2 


1.5 


1.46 


2.5 


2.45 


3.0 


2.96 


^0 


0 


0 


0 


0 


0 


0.02 


CH4 


0 


0.01 


0 


‘ 0.01 


0 


0 




0 


0.09 


0 


0.10 


0 


0.07 


c 


3.5 


3.43 


1.5 


1.48 


0 


0.06 


OH 


0 


0 


0 


0.01 


0 


0.01 





Comparison of product composition, based upon arbitrar7 and equilibrium 
relations. 



- A6 - 



RLSULTS OF THE CALCULATIONS 

Table III is an example of the results obtained by Brinkley 
and Wilson (31) using the method of approach discussed in the previous 
section. The initial and final densities, and are in grams/cubic 

centimeter. The detonation velocity, D, and the mass velocity, u^^, are 
in meters per second. The temperature is in degrees Kelvin and the 
pressure in kilobars (10^ dynes/cm^). The composition is in gram moles 
of the individual constituent per gram mole of explosive. 

A tabulation of the results obtained by Jones (27) using TNT 
- 1.5 is included for comparison. Jones' final composition gives 
larger values for CO 2 and CH^ at the expense of H 2 O and CO. The avail- 
able experimental data on detonation product composition indicate much 
smaller quantities of methane, however, it is entirely possible that 
formation may have occurred during the cooling process. The equilibrium 
of reactions involving CO and CO 2 are also sensitive to temperature changes. 
On the whole, it is reasonable to suppose that the compositions given 
by the Brinkley-Wilson procedure are a closer approach to the actual 
conditions. This is also suggested by the better agreement with the 
observed (32) detonation velocity of 6620 meter/sec. for ein initial 
density of I.5, 

The right-hand column contains results based on the arbitrary 
decomposition equations. It appears that the detonation velocity, vdiich 
is the only property that can be measured accurately, is comparatively 
insensitive to the composition. Although the three methods give results 
of the same order of magnitude, it is felt that the equilibrium approach 
of Brinkley and Wilson is preferable. 
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A comparison of observed and calculated detonation velocities 
of some common explosives is shown in Figures 7» 8 and 9. The observed 
velocities are those reported by the Explosives Research Laboratory (32) > 
Bruceton, Pa. The calculated values are based on the method of Brinkley 
and Wilson with consideration of the seven primary equilibria. The 
overall agreement is considered to be very good. 

Deviations may be due to errors in the obsenred quantities, 
however, these are usually accurate to within 3 ~ 4$, or to errors in 
the heat of formation of the explosive or to inadequacies of the method 
due to lack of accurate data in regard to covolume factors and equilibrium 
constants under the extreme conditions. 

The heats of formation of the e^qplosives considered in Figures 
7, 8 and 9 are based on the values reported by Schmidt (33) in 1924. 
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Table III 



Calculated Values of the Detonation Velocity 
and of the Properties of the Detonation Wave 
TNT 



ft 


1.0 


1.1 


1.2 


1.3 


1.4 


1.5 


1.6 


* 

1.5 


1.6 


Pi 


1.347 


1.464 


1.578 


I069O 


1.801 


1.911 


2.019 


1.94 




D 


4870 


5190 


5540 


5920 


6340 


6820 


7290 


7720 


7210 


1^1 


1250 


1290 


1330 


1370 


lUO 


1460 


1510 


1750 




Tl 


3260 


3270 


3270 


3240 


3210 


3190 


3170 


3460 


2880 


Pi 


61.0 


73.6 


88.6 


105.0 


125.4 


151.6 


177.8 


205 


173 


Composition 
H2 .09 


.05 


.02 


0 

0 

K* 


0 


0 


0 


0 


0 


CO2 


.18 


.21 


.25 


.31 


.38 


.46 


.54 


2.64 


0 


CO 


3.38 


3.22 


3.08 


2.92 


2.75 


2.58 


2.42 


.06 


3.5 


H2O 


2.26 


2.36 


2.42 


2.46 


2.48 


2.48 


2.48 


1.10 


2.5 


N2 


1.46 


1.47 


1.48 


1.49 


1.49 


1.48 


1.48 


1.50 


1.5 


NO 


0 


0 


0 


0 


.01 


.02 


.02 


0 


0 


CH4 


.01 


0 


0 


0 


0 


0 


0 


.70 


0 


NH3 


.09 


.06 


.04 


.02 


.01 


.01 


.01 


0 


0 


C 


3.43 


3.57 


3.67 


2.77 


3.87 


3.96 


4.04 


3.84 


3.5 



* Based on Jones* equation of state. 

^ Based on arbitrary decomposition equations. 
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APPiiJjDIX A 



RELA'llON OF DETONATION VELOCITY TO c 4- u AT VARIOUS POINTS 
ON THE RAMKINE-HUQQNIQT CURVE 



In a previous section, it was shown that at the point B, 
Figure 3> the detonation velocity was equal to c u in accordance with 
the Chapman-Jouguet condition. It is of interest to determine the 

relation at points other than B. Figure 3 is reproduced below for 

■ 

convenient reference. 




Figure 3 

Generalized Rankine-Hugoniot Curve 



It is well established that the sonic velocity is related to 
the slope of an isentropic PV curve by 

oi 



.2 = V? (- 



From Equations 4. and 5> we have the relation 

*0 *1 

PqVq being evaluated at the point A and Pi^V]^ on the R-H curve. From 

2 

these two equations, it appears that the difference between and 







% 

'H; ^ 







- 53 - 



2 

(D-Uj_) at various points on the R-H curve is 

L 2 I 1 ] - 

I 3Vi/3 



ol - (D-u^)^ = vf 



V -V 



( 1 ) 



The slopes of the isentrope and R-H curves at any point (P, V) are 
related by the equation 



-[(H). (^)J 



R 

where the subscript R indicates differentiation along the R-H curve, 
Rearranging and introducing the Maxwell relation | 
we have 



1 

II 

i 




\ 3S/y 


\ dvi 



I av/s \ 



R 



av/p I avjp 



(2) 



From the energy Equation 6 for points on the R-H curve, we have 



(HI 



+Po) 



Substituting for { ) its equivalent given by the first and second 

\ 3^1/r 



laws of thermodynamics yields 





. ^0-Vl 


"Pl-Po ( 9Pi\~ 


1 


1 


_ v'l 'r 



(3) 



The sign of quantity on the left-hand side of this equation depends upon 
the relative slopes of the line joining the Initial and final states and 
the tangent to the curve at the point representing the final state. Two 
extremes are shown in the figure below, bote that at point B the two 



lines coincide 




f- 
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P 




R-H curve 



V 

Figure IQ 

Slope Relations on the R-H Curves, 



If the R-H curve is concave 




>0, the difference in slopes is 



negative above B and positive below B. It is certainly reasonable to 
e;q)ect that the R-H curve should have this type curvature for the 
majority of normal fluids. Solving Equation 3 lor the quantity within 
the square brackets and subtracting from Equation 2, evaluated at the 
point permits the velocity difference written in Equation 1 to 

be ejqjressed in terms of variables of known sign. That is 




2 

The terms within the square brackets are always positive as is 
therefore the right-hand side is positive for points on the R-H curve 
above B (where ( dS^/ is negative) and negative for points below 

B, that is where ( is positive. Hence, as was to be shown, 

D is less than when the final state (PiVi) is above B and is greater 

than dQ_4“U2_ when the final state is below B, 
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APPENDIX B 
SECONDARY EFFECTS 

The theoretical treatment of the detonation of solid e^loslves 
was most conveniently accomplished by neglecting certain effects such as 
charge diameter and unstable detonation velocities that must be considered 
In any practical application of the theory. Fundamental investigations 
of these effects and their ramifications have been made and severed 
theories have been proposed in an effort to explain the observed 
phenomena. Since a detailed description of the proposed theories sheds 
little light on the fxindamentals of the detonation process. It is believed 
advisable, for the purpose at hand, to regard them as "second-ortler” 
effects and to devote to such effects a very limited treatment. 

The effect of charge diameter on the detonation velocity has 
been observed since the earliest ejqieriments were conducted. This effect 
is observed as a decrease in the detonation velocity when an explosive is 
confined in a tube of small diauneter. Fortunately, the variation ceases 
iidien a certain reasonable diameter is reached, thus the investigation of 
other parameters can be carried out without the introduction of an 
additional variable. By the same token, few investigators have under- 
taken detailed studies of its cause and until recently a generally 
accepted e^qplanation was not available, A theory proposed by Jones (34) 
accounts for the observed variation on the basis of lateral losses. 

During the finite reaction time, the high pressure with the reaction zone 
is assumed to cause an appreciable lateral expansion with a resultant 
drop in pressure, temperature and detonation velocity. In other wjrds, 
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not all of the energy of the e^qplosive is available to maintain the ideal 
detonation velocity due to lateral dissipation. It can be shown by a 
detailed mathematical analysis that the effect becomes significant when 
the radius of the explosive charge is of the same order of magnitude as 
the reaction zone. 

Another phenomenon observed frequently in the detonation of 
solid explosives is the transition from an unstable detonation velocity 
(one that is either too high or too low) to the stable value. Whenever 
an explosive charge is initiated by a priming charge having a different 
detonation rate, this transition must occur. A theory of the building-up 
rate as a function of the reaction zone length has been proposed (7) and 
leads to the conclusion that the ideal velocity is approached asymtotically, 
half the difference between the initial and final value being covered 
in a few reaction zone lengths. 

If a detonation wave experiences such severe losses that the 
wave can no longer propagate at the ideal velocity, the phenanenon 
referred to as "fading” is observed. As in the case of transition from 
one velocity to another, the length of the reaction zone is a critical 
feature and theoretical treatments of the phenomenon are based upon the 
cumulative effects i*esulting from increased reaction zone lengths. 

In summary, it might be obsez^red that if one possessed sufficient 
infomation concerning the chemical reaction kinetics to predict the 
length of the reaction zone under any conditions of pressure, temperature, 
etc., a more rational approach to many of the technical problems con- 
fronting the explosives engineer could be made. 
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